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Summary 

An  accurate  estimate  of  bit-error  rate  (BER)  in  communication  systems  is  essential  to  maintain  ac¬ 
curate  operation  of  the  system.  Currently  BER  estimate  is  assessed  through  off-line  measurements 
which  are  typically  carried  out  using  bulky  and  expensive  equipment.  In  addition  to  being  a  time 
consuming  measurement  method,  it  requires  interruption  of  traffic,  a  problem  not  allowed  in  military 
communication  links.  To  tackle  this  problem,  we  have  implemented  a  test-bed  experimental  system 
for  real  time  BER  estimator  using  an  artificial  Neural  Network  (NN)  approach  which  was  trained  off¬ 
line,  programmed  in  C  code  and  implemented  using  a  PIC32  microcontroller.  To  create  pseudo  error 
conditions,  without  transmission  of  any  additional  data,  an  additional  pseudo  error  optical  path  (at 
receiver)  is  created  splitting  the  laser  link  into  two  separate  paths.  Measurement  results  of  BER  are  in 
relatively  close  agreement  with  results  obtained  with  dedicated  BER  equipment  and  results  found  in 
literature.  The  idea  of  on-line  BER  estimate  can  provide  feedback  for  an  adaptive  control  mechanism. 
For  example,  averaged  BER  signal  (analog  value)  can  be  used  for  a  pointing  telescope  control. 
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1  Introduction 


The  goal  of  this  project  has  been  to  investigate  the  possibility  of  using  on-line  Bit  Error  Rate  (BER) 
measurements  in  laser  communication  links.  The  main  reasons  for  the  investigation  are: 

•  Existing  BER  measuring  equipment  is  bulky  and  expensive. 

•  To  measure  BER  with  this  equipment,  communication  traffic  has  to  be  interrupted. 

We  proposed  to  investigate  the  possibility  of  on-line  BER  measurement  using  an  artificial  Neural 
Network  (NN)  approach,  where  the  BER  unit  estimator  is  an  integral  part  of  the  receiver.  As  a 
result  of  this  investigation,  a  BER  monitoring  unit  was  built  using  a  pseudo  error  approach.  Figure  1 
presents  a  simplified  block  diagram  of  the  experimental  system.  The  transmitter  consists  of  a  digital 
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Figure  1 :  Block  diagram  of  experimental  BER  monitor  system. 

pattern  generator  with  scrambler  circuit,  a  1.5/im  laser  source  with  a  Mach-Zender  modulator,  and 
a  launching  optics  system,  composed  of  collimating  lenses  and  relay  mirrors.  The  receiver  consists 
of  focusing  optics  elements,  a  50/50  beam  splitter,  two  photodetectors,  a  bit  synchronizer  circuit,  a 
descrambler  and  a  BER  unit.  A  photograph  of  the  experimental  set  up  is  shown  in  Fig.  2.  This  report 
describes  the  two  essential  parts  of  the  system:  electrical  and  optical.  The  first  part  of  the  report 
concerns  the  electrical  part,  and  the  second,  the  optical  part.  In  part  three,  some  recommendation  and 
concluding  remarks  for  further  work  are  given,  and  part  four  (appendix)  lists  the  program  code  and 
measurement  table. 
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Figure  2:  Experimental  system  set  up. 


2  Electronics  Part 

2.1  The  Electronics  Transmitter 

The  transmitting  part  of  the  experimental  set  up  consists  of  pattern  generator  (HP  7084 IB),  an 
interfacing  circuitry,  a  10  bit  scrambler,  and  Mach-Zender  Modulator  (JDSU  Uniphase  Part  No. 
10020463-100).  Figure  3  shows  a  simplified  block  diagram  of  the  transmitter  together  with  the  trans¬ 
mitting  optics. 


Figure  3:  Simplified  block  diagram  of  transmitter. 


2.2  The  Electronics  Receiver 

A  simplified  block  diagram  of  the  receiver  is  shown  in  Fig.  4.  It  consists  of  an  optimal  and  pseudo¬ 
receiver,  a  synchronizer,  and  a  comparator.  During  propagation,  the  optical  signal  is  exposed  to 
different  types  of  disturbances  and  distortion.  To  estimate  the  error  of  the  received  signal,  the  optical 
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link  is  divided  into  two  paths:  a  main  link  A  (no  additional  attenuation)  and  pseudo  path  (with  some 
additional  attenuation  at  the  input  of  the  receiver).  The  idea  of  using  the  parallel  pseudo  path  is 
to  estimate  the  BER  in  the  main  path,  before  it  goes  out  of  synchronization.  For  example,  if  the 
introduced  attenuation  of  the  pseudo  path  is  two  times  higher  than  the  attenuation  of  the  main  path, 
then  the  signal,  in  the  pseudo  path,  is  attenuated  an  extra  3  dB.  This  means  that  the  pseudo  path  of  the 
link  will  lose  synchronization  before  the  main  path.  Comparing  synchronized  data  of  the  main  link 
with  a  non-synchronized  and  additionally  attenuated  data  of  pseudo  link,  we  can  get  an  estimate  of 
the  BER  in  main  optical  link.  With  good  statistics  of  the  propagation  conditions,  an  artificial  Neural 
Network  (NN)  can  be  used  as  a  BER  estimator. 


Optimal 

Receiver 


Pseudo 

Receiver 


Bit 

Synch 


Comparator 


BER 


Figure  4:  Simplified  block  diagram  of  receiver. 


2.3  The  Pseudo-Error  Detector  Circuit 

To  verify  the  pseudo  error  BER  estimate  idea,  we  have  built  the  experimental  set  up  shown  in  Fig.  5. 
Here,  the  main  receiving  path  (A)  consists  of  an  optimal  input  filter,  an  attenuator  Al,  a  threshold 
detector,  an  interface  circuit,  a  bit  synchronizer,  a  descrambler  and  the  HP70842  B  error  performance 
analyzer.  The  pseudo  receiving  path  (B)  consists  of  a  non-optimal  input  filter,  an  attenuator  A2, 
a  threshold  detector,  an  error  detector  circuitry,  and  a  PIC  microcontroller  with  the  NN.  The  error 
detector  circuit  consists  of  a  delay  line  (for  fine  tuning),  a  digital  comparator  (XOR  gate),  a  lowT  pass 
filter,  and  an  analog  multiplier  with  amplifier. 


Figure  5:  Block  diagram  of  implemented  system  (Receiver). 

The  main  path  (A)  is  synchronized  and  evaluated  in  the  HP  BER  analyzer.  The  BER  data  are 
recorded  and  evaluated  for  different  levels  of  attenuation  of  the  optical  link.  Synchronized  data  of 
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A  link  are  compared  with  non-synchronized  data  of  pseudo  link  B.  In  the  case  of  perfect  synchro¬ 
nization,  A  =  B,  the  output  of  the  digital  comparator  (XOR  gate)  is  logical  zero  (thus  BER  =  0).  As 
soon  as  the  ideal  propagation  conditions  change,  there  is  a  higher  probability  that  the  pseudo  path 
will  show  more  errors  (because  of  additionally  introduced  attenuation),  and  the  output  of  the  XOR 
detector  will  indicate  error  pulses.  For  statistical  proposes,  this  signal  is  averaged  with  a  low  pass 
filter  and  squared  to  get  the  variance.  Both,  Xi  and  X2  signals  are  fed  as  input  to  the  programmable 
NN  embedded  in  the  PIC  microcontroller.  The  complete  schematic  of  the  implemented  error  detector 
circuit  is  shown  in  Fig.  6. 
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Figure  6:  Schematic  of  a  pseudo-error  detector  circuit. 


2.4  The  Microcontroller 

Our  previous  experiments  used  Motorolas  HCS12  microcontroller.  Working  with  this  controller  was 
tedious  and  required  a  significant  amount  of  time  to  become  familiar  with  it.  Our  decision  to  switch 
to  a  PIC32  microcontroller  was  based  on  the  following  facts  [1]: 

•  The  PIC32  microcontroller  kit  contains  everything  needed  to  develop,  program,  debug  and  run 
codes  as  shown  in  a  the  block  diagram  of  Figure  7. 

•  It  is  easy  to  interface  with  an  USB  cable. 

•  It  is  easy  to  learn  how  to  program. 

•  It  is  low  cost. 
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Figure  7:  Block  diagram  of  PIC32  microcontroller. 


A  block  diagram  of  the  PIC32  microcontroller  is  shown  in  Figure  7.  The  first  step  in  the  software 
implementation  of  the  NN,  using  P1C32,  was  collection  and  evaluation  of  reliable  BER  data,  taken 
with  an  HP  BER  meter  setup  (HP70841B  and  HP70842B). 

In  the  time  period  of  two  weeks,  we  made  a  number  of  measurements,  at  room  temperature,  for 
various  levels  of  attenuation  of  an  optical  link,  in  order  to  achieve  good  statistics.  These  data  are 
averaged  and  classified  in  groups  according  to  the  level  of  a  BER  measured  with  the  HP  BER  meter. 
We  accepted  1 0-3  to  be  the  highest  acceptable  level  of  BER,  and  1 0-8  to  be  the  lowest.  The  list  of  five 
classified  groups  of  BER  measurement  is  shown  in  Appendix  A  of  this  report.  These  measurements 
were  taken  at  different  levels  of  attenuation  of  the  optical  signal.  Figure  8  shows  the  output  of  the 
error  detector  (XOR  gate)  when  the  level  of  the  error  is  1 0-8. 


Figure  8:  Output  of  the  error  detector  in  the  presence  of  attenuation  of  an  optical  signal.  Signal  1 
input  1  in  X-OR  gate  (see  figure  5),  signal  2  input  2  in  X-OR  gate.  Signal  3  output  of  X-OR  gate, 
signal  4  X2,  DC  input  in  controller,  signal  after  RC  filter 


2.5  Software  Implementation  of  NN 

For  our  proposes,  we  chose  a  two-layer  perceptron.  In  addition  to  the  two  inputs  signals  (Xj  and  X2 
in  Fig.  6),  the  perceptron  has  7  hidden  nodes  (layer  #1 ),  and  one  output  (layer  #  2).  Both  layers  have 
programming  capabilities  of  weights  and  biases.  Figure  9  shows  the  block  diagram  of  the  perceptron- 
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based  NN  used  in  our  application,  and  Figure  10  shows  a  block  diagram  of  a  single  neuron.  Its  output 
is  calculated  as  y  =  F(wx  +  b). 


Figure  9:  Block  diagram  of  adaptive  perceptron  (Adeline). 

The  Transfer  function  of  the  perceptron  adopted  is  of  the  form  y  =  1/(1  +  e_n),  and  its  graph  is 
shown  in  Figure  1 1 . 


Figure  10:  Single  input  neuron  as  a  building  block  of  the  adaptive  perceptron  in  Figure  9.  X  Input, 
W  Weight,  b  Bias,  F  -  Transfer  function,  Y  Output. 

Figure  9  consists  of  three  blocks:  a  feed  forward  perceptron,  a  training  block,  and  a  training 
algorithm.  The  procedure  to  train  the  is  as  follows:  first,  it  is  necessary  to  form  a  training  set.  This 
set  is  formed  in  such  a  way  that  the  values  of  input  signals  X!  and  X2  are  entered  (stored)  in  memory 
location  P.  At  the  same  time  the  BER  is  measured  with  an  HP  BER  meter  and  this  reading  is  recorded 
and  entered  into  memory  location  T.  After  the  training  set  is  formed,  the  values  of  the  memorized 
signals  Xi  and  X2  are  fed  back  to  the  inputs  of  the  NN.  For  these  values  of  input  signals,  the  output 
of  the  NN,  Y21,  is  subtracted  from  the  training  signal  t.  If  k2i  =  t ,  then  the  error  signal  Error 
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=  0.  This  means  that  the  coefficients  and  biases  on  the  NN  are  already  adjusted.  In  the  case  that 
error  signals  exist  (positive  or  negative),  the  training  block  for  adjustment  of  coefficients  and  biases 
generates  signals  wnew  and  bnew  (for  coefficients  adjustment,  see  Figure  9).  For  every  input  signal, 
this  procedure  is  repeated.  In  our  case,  measurements  were  made  for  six  different  levels  of  input 
signals  Xi  and  X2  (in  our  experiment,  maximum  and  minimum  acceptable  values  of  BER  are  10-3 
and  10-8  respectively).  In  Appendix  A,  the  measured  values  are  given  of  a  BER  at  different  levels 


Figure  1 1 :  Sigmoid  transfer  function. 

of  input  signals.  Their  mean  values  and  respective  variances  are  used  as  training  data  for  the  NN.  A 
weight  correction  of  layer  coefficients  is  based  on  an  algorithm  proposed  in  reference  [4],  and  its  C 
code  is  enclosed  in  Appendix  B. 

2.6  Measurement  Results 

The  system  curve.  BER  as  a  function  of  attenuation  level  of  the  optical  signal  is  shown  in  Figure  12. 
As  can  be  seen,  it  is  very  steep  and  it  is  in  general  agreement  with  the  result  published  in  the  litera¬ 
ture  [2]. 

All  the  most  relevant  equipment  purchased  for  the  electronics  system  are  summarized  as  follows: 

•  2.5GHz  High  Speed  Oscilloscope:  in  order  to  monitor  the  quality  of  the  signal,  a  high  perfor¬ 
mance,  high  speed,  digital  scope  from  Agilent  (Infiniium  Series  No.  DSO90254A)  was  used.  It 
is  crucial  for  the  system  to  work  properly  that  the  received  signal  displays  sharp  edges,  i.e  tran¬ 
sitions  between  the  low  values  (logical  ’0’s)  and  high  values  (logical  ’l’s).  That  requires  fast 
response  equipment.  For  instance,  to  detect  a  signal  with  a  rise  time  of  1  nsec,  it  would  roughly 
require  a  system  with  a  1GHz  bandwidth  capability.  This  scope  provided  2.5GHz  bandwidth 
with  20Gsamples/sec  sampling  speed,  which  was  more  than  adequate. 

•  EZJIT  jitter  Analysis  Software:  this  software  (DS090000A-004)  was  part  of  the  extended 
package  purchased  together  with  the  Infiniium  scope,  and  it  is  an  essential  tool  for  jitter  signal 
analysis. 

•  High  Performance  Probe:  a  high  performance,  differential  browse  probe  (E2675A),  from  Ag¬ 
ilent,  with  a  high  impedance  adapter  (E2697A),  and  optionals  (1 131A,E2669A),  was  purchased 
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Figure  12:  (a)  Results  of  measurement,  and  (b)  Results  from  reference  [2], 


to  assist  with  signal  monitoring  at  points  of  the  system  where  any  disturbance  of  the  node  would 
be  critical.  One  such  a  probe  was  necessary,  for  its  capacity  of  “tapping”  out  a  sample  of  the 
the  signal  without  changing  the  characteristics  of  the  information  flow  at  that  point  on. 

•  High  Precision  Delay  Line:  to  compensate  for  signal  “de-phasing”,  one  of  the  most  fundamen¬ 
tal  requirements  in  the  lab  was  a  precision  delay  line.  For  that,  a  programmable  delay  line  from 
Colby  Instruments  (HPDL-100A-10.23NS)  was  purchased.  This  delay  line  has  the  capacity  to 
add  delays  in  as  little  as  1  Opsec  increments  within  a  1  Onsec  range. 

•  Mixed  Signal  Scope:  on  those  points  of  the  circuit  where  speed  was  not  an  issue,  and  which 
required  constant  monitoring,  a  low  speed,  low  bandwidth,  scope  (MSO7034A)  and  probes 
(10074C),  from  Agilent,  was  used.  This  is  a  much  rudimentary  equipment  compared  to  the 
(DSO90254A)  scope,  but  nonetheless  essential. 

•  Free  Space  Broadband  Detector:  several  photodetectors  were  used  in  the  project.  For  those 
specific  measurements  where  bandwidth  was  a  requirement,  i.e  sharp  edge  monitoring,  a  high 
bandwidth  detector  was  necessary  to  match,  in  fact,  surpass  the  bandwidth  of  the  fast  scope 
(DSO90254A).  For  that,  a  25GHz  broadband  photodetector  from  NewFocus  (1437)  was  pur¬ 
chased. 

•  Power  Supplies:  several  high  quality  power  supplies  were  necessary  to  bias  all  the  electronic 
circuitry,  at  several  voltage  and  power  levels.  The  best  choice  was  to  use  power  supplies  from 
Agilent  (E3616A,E3620A,E3610A). 

•  Buffers  and  Voltage  Translators:  to  assist  with  signal  distribution,  regeneration,  and  volt¬ 
age  level  translation,  required  for  interfacing  electronic  circuits  of  different  voltage  standards, 
several  50ohm  line  drivers  (fannout  buffers)  and  translators,  from  Pulse  Research  Lab  (PRL- 
420ND,PRL-450ND,  PRL-414B)  were  purchased. 
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3  BER  Estimate  for  Adaptive  Optics  Control 

Laser  communication  systems,  such  as  in  airbom  and  submarine  platforms,  use  line  of  sight  optical 
links.  In  such  links  it  is  important  that  communication  link  remains  invisible.  Link  between  platforms 
must  be  reliable  in  the  presence  of  a  variety  of  disturbances  such  as  atmospheric  and  hidrospheric 
perturbations,  movements,  bums,  etc.  To  overcome  this  obstacles,  a  reliable  control  of  communication 
link  must  be  achieved.  In  our  opinion,  BER  estimate  of  the  laser  link  can  be  used  for  this  goal. 
Figure  13  presents  a  simple  block  diagram  of  the  proposed  control  system  based  on  BER  estimate. 


BER 

Figure  13:  Block  diagram  of  proposed  control  system. 

A  BER  estimate  is  detected  in  point  3  (the  output  of  the  XOR  gate  in  Figure  5)  and  is  modulated 
by  a  laser  modulator,  then  reflected  by  a  mirror  and  detected  by  four  quadrant  diodes.  Quadrant 
diode  detectors  consists  of  4  photo  sensitive  diodes  and  the  outputs  produce  voltages  proportional  to 
the  amount  of  photon  flux.  When  a  laser  beam  is  focused  in  the  middle  of  the  detector,  then  every 
diode  has  one  quarter  of  its  light,  and  the  voltage  levels  at  the  outputs  of  the  diodes  are  equal.  As  the 
beam  moves  off-center,  the  proportion  of  light  on  each  diode  changes.  As  a  consequence,  their  output 
potentials  change  accordingly.  These  changes  in  voltages  are  detected  by  BER  electronics  (digital 
BER  signal  is  averaged  to  obtain  a  nearly  DC  level  signal)  and  processed  in  micro-controller.  The 
differences  in  the  diodes  potentials  determine  the  magnitude  and  location  of  error  in  the  beam  spot 
position.  Beam  steering  can  be  accomplished  with  a  mirror  mounted  on  2  DC  motors. 


4  The  Optical  Link 

The  optical  link  is  the  part  of  the  system  where  the  information  was  carried  optically  using  a  laser 
beam.  It  relates  to  the  generation,  modulation,  propagation  and  detection  of  the  optical  signal.  During 
the  project,  several  changes  in  the  optics  were  made  to  enhance  the  performance  of  the  link.  The 
following  sections  summarize  the  evolution  of  the  system  from  one  year  to  the  next. 
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In  the  first  year,  the  most  basic  optical  link  was  assembled.  It  was  a  single  channel  system,  where 
the  optical  signal  was  propagated  all  the  way  up  to  the  roof  of  the  CHTM  building  facility  and  back. 
Figure  14  shows  the  schematics  of  the  link. 


Opto-electronics  System 


RF  IN 


Figure  14:  Schematics  of  the  basic  optical  link  system  implemented  in  FY#1:  opto-electronics, 
launching/receiving  optics  and  receiver  system.  Blue  lines  indicate  optical  fiber  paths.  Red  lines 
indicate  light  in  air. 

The  system  is  divided  in  three  separate  sections:  (1)  the  opto-electronics  system,  where  the  op¬ 
tical  signal  is  generated,  modulated,  and  amplified,  (2)  the  launching  system,  where  the  amplified 
modulated  light  is  propagated  in  air  and  directed,  through  lenses  and  relay  mirrors  to  the  third  sec¬ 
tion,  and  (3)  the  receiver  system,  where  it  is  collected  and  imagined  onto  a  photodetector  (PD)  to  be 
transformed  back  into  an  RF  modulated  signal. 

4.1  The  Opto-electronic  System 

In  the  opto-electronics  system,  a  Laser  Diode  (LD)  light  source  is  directly  connected  to  an  optical 
Amplitude  Modulator  (AM),  which  is  interfaced  to  the  electronics  of  the  system  via  an  RF  input 
signal  port.  The  output  of  the  modulator  goes  to  an  Erbium  Doped  Fiber  Amplifier  (EDFA),  for 
amplification  using  standard  Single  Mode  Fibers  (SMF).  The  most  important  components  purchased 
for  this  system  are  listed  below: 

•  The  Laser  Diode  (LD)  Source:  a  PM  fiber  coupled,  20mW,  1550nm  center  wavelength,  CW 
DFB  laser  diode,  from  JDSU  Uniphase  (CQF935/508  Series).  The  LD  was  driven  using  a 
Spectra  Diode  Labs  laser  diode  driver  (SDL-800),  and  operated  at  an  average  output  power  of 
15mW. 


•  The  Amplitude  Modulator  (AM):  a  lOGb/sec,  “bias-ready”  intensity  modulator  from  JDSU 
Uniphase.  The  modulator  was  coupled  with  PM  fiber  at  the  input,  and  directly  attached  to  the 
laser  source  via  an  FC/UPC  mating  sleeve  connector.  The  output  fiber  was  standard  SMF. 

•  The  Optical  Amplifier:  an  Erbium-Doped  Fiber  Amplifier  (EDFA)  from  IPG  Photonics  (Model 
EAR-5K-C).  The  amplifier  operates  with  a  minimum  input  power  of -5dBm,  with  was  suitable 
to  the  low  power  output  level  coming  out  of  the  modulator. 

4.2  The  Launching  System 

In  the  launching  system,  or  emitter,  the  light  coming  from  the  EDFA  is  collimated  using  a  fiber  ready, 
FC  female  type  collimator  (COL),  followed  by  a  beam  expander.  From  this  point  on,  a  series  of  metal 
film,  flat  mirrors  (LM  i,  LM2)  relayed  the  5cm  diameter  optical  beam  through  the  buildings  roof,  back 
into  the  optical  table  (REF,  RMi,  RM2,  RM*),  where  a  single  biconvex  lens  (1Q)  focused  the  light 
into  an  InGaAs  type  photodetector  (PD).  In  the  second  year,  an  improved  system  for  the  emitter  was 
implemented.  In  comparison  to  the  system  used  in  the  previous  year,  this  new  system  offered  the 
following  advantages: 

•  Splitting  of  the  emitted  optical  beam:  as  mentioned  in  the  original  proposal,  splitting  the  laser 
into  multiple  beams  enhances  the  performance  of  the  system,  as  it  helps  minimize  the  effects 
of  scintillation  [3,  5,  6],  which  results  from  propagation  in  the  atmosphere.  In  this  new  system, 
the  beam  was  split  into  4  separate  channels  using  a  high  power  fiber  optic  coupler. 

•  Attachment  of  the  emitter  to  the  beam  director:  an  important  innovation  of  the  new  optical 
system  is  the  assembly  of  the  whole  launching  optics  on  the  frame  of  the  beam  director  located 
on  the  buildings  roof.  In  this  way,  the  entire  launching  system  was  attached  to  the  frame  of 
rotation  of  the  beam  director,  and  alignment  of  the  emitter  was  no  longer  necessary  every  time 
the  beam  director  was  rotated. 

•  Use  of  larger  aperture  launching  mirrors:  Compared  to  the  3  diameter  mirrors  used  in  the 
previous  year,  new  6  diameter  flat  mirrors  allowed  the  possibility  of  placing  the  target  at  much 
farther  distances,  as  it  can  be  used  to  increase  the  Rayleigh  range  of  the  propagating  beam  . 

Schematics  of  the  new  emitter  system  are  shown  in  Figure  15.  It  starts  with  the  amplified  optical  sig¬ 
nal,  which  comes  out  of  the  Erbium-doped  fiber  amplifier.  The  collimated  beam  is  coupled  to  a  section 
of  Single-Mode  Fiber  (SMF)  using  an  XYZ  fiber  launching  system  (Thorlabs,  Model  No.MBT612), 
as  seen  in  Figure  16.  The  other  end  of  the  fiber  is  connected  to  a  high  power  1  x  4  optical  coupler 
(NEPTEC  No.  CC-FA14BO0419E)  that  splits  out  the  beam  into  4  separate  channels,  where  the  light 
propagates  along  4  separate  sections  of  SM  fibers.  The  output  fibers  are  taken  directly  to  the  roof, 
inside  the  beam  director  (Fig.  20).  There,  they  are  mounted  on  90  degree  racks,  together  with  15cm 
focal  length  lenses  to  form  a  collimating  assembly  (see  inlet  in  Fig.  17).  The  racks  are  placed  right 
below  4  6“  metal  flat  mirrors  which  are  then  used  to  relay  the  beams  to  the  reflecting  target  (see  REF 
in  Fig.  14). 
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Figure  16:  Fiber  launching  system,  where  a  collimated  amplified  beam  (left)  is  coupled  into  a  SM 
fiber  using  a  microscope  objective  mounted  on  a  3  axis  positioning  stage  (right). 


All  the  most  relevant  parts  purchased  for  this  system  are  summarized  as  follows: 

•  Launching  Mirrors  (LM^):  4  6“dia.  fused  silica,  protected  silver,  flat  mirrors  (CVI  Custom 
Specified)  were  at  the  final  stage  of  the  emitter,  launching  the  collimated  beams  to  the  reflecting 
mirror  (REF)  which  pointed  them  toward  the  receiver’s  mirror. 

•  Beam  Splitter  (BS):  a  1x4  high  power  beam  splitter,  coupled  to  the  output  of  the  EDFA,  split 
the  beam  into  4  separate,  equal  power  channels.  The  beam  splitter  used  was  from  NEPTEC 
Optical  Solutions  (P/N  CC-SA14B0Q419). 
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•  Optical  Fiber  Patchcords:  4  8  meters  long  optical  fiber  patchcords  carried  the  signal  all  the 
way  up  to  the  building’s  roof,  where  the  beam-director  was  placed.  The  patchcords  were  SMF 
standard,  connected  to  the  ports  of  the  BS  via  FS/APC  mating  sleeve  connectors.  The  APC 
types  of  connectors  come  with  an  angle  at  the  fiber  end  to  minimize  reflections  back  into  the 
system,  which  would  be  damaging  at  those  power  levels. 

•  Collimators  (CLi_4):  four  2“  dia.  IR  CaF2  singlet  spherical  plano-convex  lenses  (CVI  PLCX- 
50.8-83.8-CFIR-1550)  placed  at  the  collimation  stages  directly  above  the  fiber  ends,  inside  the 
beam  director,  provided  collimation  to  the  output  beams  in  air. 

4.3  The  Receiver  System 

In  the  third,  and  final  year  of  the  project,  an  improvement  was  made  on  the  Receiver  system.  Figure  1 8 
shows  the  schematics  of  the  system.  In  comparison  to  the  system  of  the  previous  year,  this  new  system 
provided: 

•  Addition  of  a  second  channel  to  the  received  signal:  using  a  polarization  independent  50/50 
beam  splitter,  the  received  signal  was  split  into  two  beams,  allowing  monitoring  of  the  signals 
at  two  separate  channels,  as  required  by  the  electronics  system. 


Figure  17:  Photograph  of  the  entire  emitter  system:  the  launching  optics  inside  the  beam  director 
(background),  and  detail  of  the  fiber  collimator  assembly  (inlet). 
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•  Insertion  of  a  linear  attenuator:  an  attenuator  was  placed  before  the  beam  splitter  to  allow 
gradual  linear  attenuation  of  the  received  signal. 

•  High  performance  photodetectors:  to  equalize  the  response  of  two  channels  of  the  receiving 
system,  two  new  high  performance  photodetectors  were  placed  at  the  end  of  the  link.  These 
new  detectors,  with  bigger  active  areas,  improved  focusing,  and,  therefore,  signal  detection. 

Figure  18  shows  the  schematics  of  the  new  receiver.  It  consists  of  5  mirrors  (RMi,  and  Mi_s),  1 
attenuator,  1  50/50  beam  splitter  (BS),  two  focusing  lenses  (l0i— 2),  and  2  photodetectors  (PDi_2).  The 
first  (and  largest)  mirror  is  mounted  at  the  beam  director,  together  with  the  launching  system  mirrors 
(the  large  mirror  at  the  center  of  Fig.  17).  It  relays  the  four  incoming  beams  to  a  mirror  placed  directly 
below  it,  inside  the  lab  (Mi),  at  the  optical  bench  (same  as  RM2  in  Fig.  17).  This  mirror  relays  the 
beam  to  M2,  which,  in  turn,  relays  to  M3.  An  attenuator,  placed  right  in  front  of  M4,  allows  gradual 
linear  attenuation  of  the  optical  signal,  which  is  then  sent  to  a  50/50,  polarization  insensitive,  beam 
splitter  (BS).  At  this  point,  half  of  the  signal  is  focused  to  one  photodetector  (PDi)  via  l0i,  whereas 
the  other  half  is  focused  to  PD2  via  lo2,  after  passing  another  relay  mirror  M3. 


Figure  18:  Schematics  of  the  receiver  system.  Red  lines  indicate  light  in  air. 

The  new  purchased  parts  which  were  added  to  the  system  are  summarized  as  follows: 

•  50/  50  Beam  Splitter  (BS):  a  2“dia.  50/50  beam  splitter  from  CVI  laser  (Part.  No.  BSNP- 
1550-50-2025). 

•  Linear  Attenuator:  a  mounted  linear,  variable,  metallic  ND  filter  (Thorlabs,  Inc.  Part  No. 
NDC-100-4M). 

•  Focusing  lenses  (l0i— 2)1  2  2“  dia.,  50cm  focal  length,  biconvex  lenses  (CVI  laser,  Part.  No. 
BICX-59.8-49-3-C,  1 08- 1 0,33-20). 

•  Photodetectors  (PD!_2):  2  0.2mm-’  active  area,  wideband  amplified  InGaAs  detectors  (Thor¬ 
labs,  Inc.  Part.  No.  PDA10CF). 
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Figure  19  shows  a  TOP  VIEW  picture  of  the  actual  new  receiver  system.  The  red  lines  indicates 
the  trajectory  of  the  optical  signal  along  the  elements  of  the  receiver. 


Figure  19:  Photograph  of  the  actual  receiver  system.  Red  lines  indicate  the  light  path  in  air. 


■ 


Figure  20:  Photograph  of  the  beam  director  located  at  the  roof  of  the  CHTM  facility. 
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5  Concluding  Remarks 

Our  experimental  results  have  shown  that  real  time  BER  monitoring  system  has  potential  applications 
in  different  areas  of  digital  communication  links.  In  addition,  being  integral  part  of  the  receiver, 
signal  of  BER  estimator  can  be  used  for  control  purposes.  Our  additional  goal  was  to  implement  NN 
system  with  on-line  training  capabilities.  Unfortunately,  I,  Dr  Zrilic,  have  faced  permanent  problem 
of  keeping  qualified  graduate  students.  In  past  three  years  six  different  graduate  students  -  Alex  P. 
Braga,  Pavlov,  Behshad  Amadi,  Zhuoyao  Wang,  Kai  Du,  Mladen  Prijic  -  were  employed. 
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APPENDIX 


A  BER  Measurements 
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Table  1 :  BER  Measurements. 
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Table  2:  BER  Measurements. 
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B  C  Program 

/* 

**  main.c 
*  * 

*/ 

//  configuration  bit  settings,  Fcy=72  MHz,  Fpb=36  MHz 
#pragma  config  POSCMOD=XT,  FNOSC=PRIPLL 
//New  Oscillator  Selection  bits 

#pragma  config  FPLLIDIV=DIV_2 ,  FPLLMUL=MUL_18 ,  FPLL0DIV=DIV_1 
//Pll  input  driver,  Pll  multiplier,  Pll  output  driver 
#pragma  config  FPBDIV=DIV_2 ,  FWDTEN=OFF,  CP=OFF,  BWP=OFF 
//Peripheral  Bus  Clock  Divisor,  WatchDog  Timer  Enable  bit, 
Code-Protect  bit,  Code-Protect  bit 

#include  <p32xxxx.h> 

#include  "Explore. h" 

#include  "ADC.h" 

#include  "LCD . h" 

#include  "NeuralNetwork . h" 

main  ( ) { 

int  read,  s,  u,  v,  samp; 
float  input 0,  inputl; 
struct  neuralnetwork  x; 

initADC (AINPUTS) ;  //  initialize  the  ADC 
initLCDO;  //  initialize  the  LCD  display 

x  =  inicialization (x) ; 
x  =  inicializationWeight (x) ; 
x  =  trainNN(x); 

while (1 ) { 

inputO=ADC ( 0 ) ; 
input 1=ADC ( 1 ) ; 

x=putlnputs (x,  inputO,  inputl); 
x=propagation (x) ; 

Delayms (2000) ; 

displayFloat (x. layers [2] .output [0] ) ; 

}  //  main  loop 

} 

/* 
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**  LCDlib.c 
★  * 

*/ 

#include  <p32xxxx.h> 
#include  <plib.h> 
#include  "Explore. h" 
#include  "LCD.h" 


# define  PMDATA  PMDIN 


void  initLCD(  void) 

{ 

//  PMP  initialization 

mPMPOpen (PMP_ON  |  PMP_READ_WR I TE_EN  |  3, 

PM  P_DATA_BU  S_8  |  PMP_MODE_MASTERl  | 

PMP_WAIT_BEG_4  j  PMP_WAIT_MID_15  | 

PMP_WAIT_END_4 , 

0x0001,  //  only  PMAO  enabled 

PMP_INT_OFF) ;  //no  interrupts  used 
//  wait  for  >30  ms 
Delayms (30) ; 

//initiate  the  HD44780  display  8-bit  init  sequence 
PMPSetAddress (LCDCMD) ;  //  select  command  register 

PMPMasterWrite  (0x38)  ;  //  8-bit  int,  2  lines,  5x7 
Delayms (1);  //> 48  us 

PMPMasterWrite (0x0c) ; 

Delayms (1);  //>48  us 

PMPMasterWrite (0x01) ; 

Delayms (2);  //>1.6  ms 

PMPMasterWrite  (0x06)  ;  //  increment  cursor,  no  shift 
Delayms(2);  //>1.6  ms 
}  //  initLCD 


//ON,  no  cursor,  no  blink 
//  clear  display 


char  readLCD (  int  addr)  //  readLCD 

{ 

PMPSetAddress (addr) ;  //  select  register 
mPMPMasterReadByte  ( ) ;  //  initiate  read  sequence 
return  mPMPMasterReadByte () ;  //  read  actual  data 
} 


void  writeLCD(  int  addr,  char  c) 

{ 

while (busyLCD () ) ; 

PMPSetAddress (addr) ; 
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PMPMasterWrite (c) ; 

} 

void  put sLCD (char  *s) 

{ 

while ( *s) 
putLCD ( *s++) ; 

} 

void  displaylnt (unsigned  int  integer) 

{ 

if (integer/10000) 

putLCD (0x30  +  (integer  /  10000)); 
if ( integer/1000 ) 

putLCD(0x30  +  (integer  %  10000)  /  1000); 
if (integer/100) 

putLCD(0x30  +  ((integer  %  10000)  %  1000)  /  100); 
if (integer/10) 

putLCD ( 0x3 0  +  (((integer  %  10000)  %  1000)  %  100)  /  10  ) 
putLCD ( 0x3 0  +  (((integer  %  10000)  %  1000)  %  100)  %  10  ) 
} 

void  displayFloat (float  number) 

{ 

int  j , i ; 

float  remainder; 
j=0; 

clrLCDO  ; 

putsLCD(  "Error  ratio"  ); 
cmdLCD (  0x80  |  0x40); 
if  (number >0) 

{ 

while  (number<l) 

{ 

j-j+i; 

number=number*io ; 

} 

if  { j  >2  5) 

{ 

number=0 . 0 ; 

displaylnt ( (unsigned  int) number) ; 
putLCD (0x2e) ; 

remainder= (number  -  ((unsigned  int) number) ) *10 ; 
displaylnt ( (unsigned  int) remainder) ; 

remainder= (remainder  -  ((unsigned  int) remainder) ) *10 ; 
displaylnt ( (unsigned  int ) remainder) ; 


22 


putsLCD ( "xlO-O " ) ; 
displaylnt ( 0 ) ; 

} 

else 

{ 

displaylnt ( (unsigned  int)number) ; 
putLCD ( 0x2e) ; 

remainder= (number  -  ((unsigned  int) number) ) *10 ; 
displaylnt ( (unsigned  int) remainder) ; 

remainder= (remainder  -  ((unsigned  int) remainder) ) *10 ; 
displaylnt ( (unsigned  int) remainder) ; 
putsLCD ( "xlO-O" ) ; 
displaylnt ( j ) ; 

} 

} 

//  write  negative  number 
else 
{ 

putLCD ( 0x2 d) ; 
number= -number ; 
while  (numbercl) 

{ 

j=j+i; 

number=number*10  ; 

} 

displaylnt ( (unsigned  int)number) ; 
putLCD ( 0x2 e) ; 

remainder= (number  -  ((unsigned  int) number) ) *10 ; 
displaylnt ( (unsigned  int) remainder) ; 

remainder= (remainder  -  ((unsigned  int) remainder) ) *10 ; 
displaylnt ( (unsigned  int) remainder) ; 
putsLCD ("xlO-O") ; 
displaylnt ( j ) ; 

} 

//Delayms (2000) ; 

} 

/* 

**  ADClib.c 
★  ★ 

*/ 


#include  <p32xxxx.h> 
#include  <plib.h> 
#include  "Explore. h" 
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#include  "ADC . h 


void  initADC(int  amask) { 
ADlPCFG=amask; 
ADlCONl=OxOOEO ; 

AD1CSSL=0 ; 

AD1CON2=0 ; 

AD1CON3  =  0xlF3  F ; 
AD1CON1SET=0x8000 ; 


} 


int  readADC(int  ch) { 
ADICHSbits . CH0SA=ch ; 
ADlCONlbits . SAMP=1 ; 
while ( lADICONIbits. DONE) ; 
return  ADC1BUF0; 

} 

float  ADC  (int  port) { 
int  i  ; 

float  average,  sweep; 
sweep=0 ; 

for(i=0;  i<20;  i++) { 
sweep  +=  readADC (port) ; 
Delayms (300) ; 

} 

average=sweep/2  0 ; 
return  average ; 

} 


/  * 

**  Explore. c 


*/ 


#include  <p32xxxx.h> 

#include  <plib.h> 

#include  "explore. h" 

void  initEX16 (  void) { 

#if ndef  PIC32_STARTER_KIT 
mJTAGPortEnable (0) ; 

#endif 

SYSTEMConf igPerf ormance (  FCY) 


INTEnableSystemMultiVectoredlnt ()  ; 

LATA  =  0; 

TRISA  =  0xFF80 ; 

} 

void  _general_exception_handler (unsigned  c,  unsigned  s)  { 
while  ( 1 ) ; 

} 

/* 

**  Simple  Delay  functions 
*/ 

void  Delayms (unsigned  t)  { 

T1CON  =  0x8000; 
while  (t--) 

{ 

TMR1  =  0; 

while  (TMR1  <  FPB/1000) ; 

} 

} 

/  * 

**  NeuralNetwork . c 
*  * 

*/ 


#include 

#include 

#include 

#include 

#include 


<p32xxxx . h> 

<plib . h> 

"explore . h" 

"LCD . h" 

"NeuralNetwork . h" 


struct  neuralnetwork  inicialization ( struct  neuralnetwork  x) { 


x. layers [0] . n=2 ; 
x. layers [1] .  n=7; 
x. layers [2] . n=l ; 

x. trainSet  [0]  [0] 
x . trainSet  [0]  [0] 
x. trainSet  [0]  [1] 
x. trainSet  [1]  [0] 
x. trainSet  [1]  [0] 
x. trainSet  [1]  [1] 


[0] =73; 

[1] =36; 

[0] =0. 000000015; 
[0] =92; 

[1] =59; 

[0] =0 . 00000015 ; 
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x . trainSet 
x. trainSet 
x . trainSet 

x . trainSet 
x . trainSet 
x . trainSet 

x . trainSet 
x . trainSet 
x . trainSet 

x . trainSet 
x . trainSet 
x . trainSet 


[2] 

[0] 

[0] 

=101; 

[2] 

[0] 

[1] 

=  70; 

[2] 

[1] 

[0] 

=  0 . 0000015 

[3] 

[0] 

[0] 

=  114; 

[3] 

[0] 

[1] 

=  86; 

[3] 

[1] 

[0] 

=0 . 000015; 

[4] 

[0] 

[0] 

=  118; 

[4] 

[0] 

[1] 

=  94; 

[4] 

[1] 

[0] 

=  0 . 00015; 

[5] 

[0] 

[0] 

=  132; 

[5] 

[0] 

[1] 

=  110; 

[5] 

[1] 

[0] 

=0.0015; 

return  x; 

} 

struct  neuralnetwork  inicializationWeight ( struct  neuralnetwork  x) { 
int  s,  u,  v; 

//  initialization  weight  in  layer  0 
s=0  ; 

for(u=0;  u<x . layers [s] . n;  u++) 
for(v=0;  v<x . layers [s+1] . n- 1 ;  v++) { 
x. weight [s] [0] [v] =0.01901; 
x. weight [s] [1] [v] =0.01101; 

} 

//initialization  weight  in  layer  1 
s  =  l ; 

for(u=0;  u<x . layers [s] . n;  u++) 
for(v=0;  v<x . layers  [s  +  1]  . n;  v++) { 
x. weight [s] [u] [v] =7.201; 

} 

//initialization  bias  in  first  layer 
for(u=0;  u<x . layers [s] . n;  u++) { 
x. layers [s] .bias [u] =-4 . 2501 ; 

} 


//initialization  bias  in  second  layer 
x.layers[2] .bias [0] =-22 . 301 ; 

//  initialization  delta  weight 
for ( s=0 ;  s<NUM_LAYERS-l;  S++) 

for(v=0;  v<x . layers [s+1] . n;  v++) 
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for(u=0;  u<x. layers [s] .n;  u++)  { 

x.weightP[s]  [u]  [v]=0.0; 

} 


//  initialization  delta  bias 
for (s=l ;  s <NUM_LAYERS - 1 ;  s+  +  ) 

for(v=0;  v<x . layers  [s]  . n;  v++)  { 

x . layers [s+1] .biasp [v] =0 . 0 ; 

} 

return  x; 

} 

float  sigmoid (float  net)  { 
return  1/ (l+exp(-net) ) ; 

} 

struct  neuralnetwork  propagation ( struct  neuralnetwork  x)  { 
int  s,  v,  u; 
float  net  =  0.0; 
for (s=l ;  s<NUM_LAYERS ;  S++) 
for (v=0;  v<x . layers [s] . n ;  v++)  { 

net  =  x . layers [s]  .bias  [v]  ; 

for(u=0;  u<x. layers [s-1] .n;  u++) 
net  +=  x. layers [s-1] .output [u] *x. weight [s-1] [u] [v] ; 
x. layers [s] . output [v]  =  sigmoid (net) ; 

} 

return  x; 

} 

struct  neuralnetwork  putlnputs (struct  neuralnetwork  x, 
float  inputl,  float  input2) { 
x. layers [0] .output [0]  =  inputl; 
x . layers [ 0] . output [1]  =  input2 ; 
return  x; 

} 

struct  neuralnetwork  putlnputsTrainSet (struct 
neuralnetwork  x,  int  set) { 
int  u  ; 

for(u=0;  u<x . layers [0] .n;  u++) 

x. layers [0] . output [u]  =  x . trainSet [set] [0] [u] ; 
return  x; 

} 

struct  neuralnetwork  error  (struct  neuralnetwork  x,  int  set)  { 
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int  v,  u,  s; 
float  f,  sigmaa; 

//  error  in  output  (second)  layer 
for(  v=0;  v<x. layers [NUM_LAYERS-1] .n;  v+  +  )  { 
x. layers [NUM_LAYERS- 1] . difference [v]  = 
x. trainSet  [set]  [1]  [v]  -  x . layers [NUM_LAYERS- 1]  . output [v]  ; 

x. error  +=  (x,.  layers  [NUM_LAYERS-1]  .  difference  [v]  ) 

* (x. layers [NUM_LAYERS-1]  .difference [v]  )  ; 

x . layers [NUM_LAYERS - 1 ] . sigma [v] =x . layers [ 

NUM_LAYERS - 1 ]  .difference  [v] *x. layers [NUM_LAYERS- 1]  .output [v] * 
(1  -  x. layers [NUM_LAYERS-1] . output [v] ) ; 

} 

//  error  in  first  layer 
for (s=NUM_LAYERS-2 ;  s>=0;  s--){ 

for(u=0;  u<x . layers [s] . n;  u++)  { 

sigmaa  =  0.0; 

for(v=0;  v<x. layers [s+1] .n;  v++) { 

sigmaa  +=  x . layers [s+1] . sigma [v] *x . weight [s] [u] [v] ; 

} 

f  =  x . layers [s] . output [u] ; 

x. layers  [s]  . sigma [u]=  f * (1-f) ★sigmaa; 

} 

} 

return  x; 

} 

struct  neuralnetwork  trainNN ( struct  neuralnetwork  x) 

{ 

int  iteration,  samp,  s,  v,  u; 

f or ( iteration=0 ;  iteration<MAXITERATION;  iteration++) { 
x . error=0 ; 

for(samp=0;  samp<SET;  samp++) { 
x  =  putlnputsTrainSet (x, samp) ; 
x  =  propagation (x) ; 
x  =  error  (x,  samp) ; 
x. error  +=  0.5*x.error; 
for ( s  =  0 ;  s  <NUM_LAYERS - 1 ;  s  +  +  )  { 

for(v=0;  v<x . layers [s+1] . n;  v++)  { 

for(u=0;  u<x. layers [s] .n;  u++)  { 

x.weightP[s]  [u]  [v]  =  NI*x. layers [ 
s+1]  . sigma  [v] *x. layers [s]  . output [u] ; 

x. weight [s]  [u]  [v]  +=  x.weightP[s]  [u]  [v]  ; 

} 

x . layers  [s+1]  . biasp [v]  = 
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NI*x . layers [s+1] . sigma [v] ; 

x . layers [s+1] . bias [v]  +=  x. layers [s+1] .biasp [v] ; 

} 

} 

} 

if (x . error<MAXERROR) 
break; 

} 

return  x; 

} 
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